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Abstract: In the Color Glass Condensate formalism, charged hadron px distributions in p+p 
collisions are studied by considering an energy-dependent broadening of nucleon’s density distribu¬ 
tion. Then, in the Glasma flux tube picture, the n-particle multiplicity distributions at different 
pseudo-rapidity ranges are investigated. Both of the theoretical results show good agreement with 
the recent experimental data from ALICE and CMS at ^/s = 0.9, 2.36, 7 TeV. The predictive results 
for pt and multiplicity distributions in p-|-p and p-l-Pb collisions at the Large Hadron Collider are 
also given in this paper. 
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1. Introduction 

The measurement of charged hadron transverse momentum (px) spectrum in heavy-ion collisions 
is a very significant topic on experimental physics. Recently, the ALICE and CMS experiments at 
the Large Hadron Collider (LHC) have given the data for charged hadron px distributions in p-|-p 
collisions at ^/s = 0.9, 2.36, 7 TeV [1-2]. The experimental results show that the scaling features 
of charged hadrons previously observed at lower energies do not hold good at high energies. These 
indicate that we still can not give some basic understanding of the mechanism involved in the charged 
particle production. Thus, a systematic theoretical study on this topic is still very necessary. 

An effective theory to study charged hadron px spectrum is the Color Glass Condensate (CGC) 
approach [3]. The CGC is a state predicted by quantum chromodynamics (QCD) at high energies 
where gluons in a hadron wave function expect to statute. The cornerstone of the CGC approach is 
the existence of a hadron saturation scale Qs at which gluon recombination effects start to balance 
gluon radiation. In order to give an accurate theoretical analysis, the nucleon’s density distribution 
in position space must be considered. Thus, the saturation scale, Qs, should be considered as a 
function of the nucleon’s thickness function, Tp{b). Eurthermore, due to gluon saturation, the width 
of the gluon distribution inside a nucleon should grow with collision energy, y/s [4]. This will lead to 
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a broadening of the nucleon’s density distribution in position space as ^/s increases. Therefore, an 
energy-dependent broadening of the nucleon’s thickness function in the Qs should also be considered. 

Based on the CGC theory, many phenomenological saturation models, such as the Golec-Biernat 
and Wiisthoff (GBW) model [5], the Kharzeev, Levin and Nardi (KLN) model [6], the CGG model 
[7] and the Kovchegov, Lu and Rezaeian (KLR) model [8], are established. R should be noted that 
the KLR model is based on the Anti-de Sitter space/conformal field theory (AdS/GFT). In this 
paper, the KLN and KLR-AdS/GFT model are used because the simple analytic unintegrated gluon 
distribution functions can be obtained by a Fourier transform from the dipole amplitude of these 
two models. 

Integrating charged hadron transverse momentum distributions over px, we can obtain the im¬ 
pact parameter (b) dependent mean multiplicity, h(b). The mean multiplicity and the parameter 
controlling the size of the fluctuations, k, are two parameters that characterize negative binominal 
distribution (NBD) [9-10]. In the GGC picture, particles produced locally in the transverse plane 
are considered as correlating by approximately boost invariant flux tubes [11]. Since there are Q‘gSj_ 
such flux tubes emitting gluons, the parameter k should be proportional to Q‘sS±, where S± is the 
transverse overlap area of two collision hadrons. This indicates that the parameter k in the NBD is 
also a function of the thickness function, Tp{b), and influenced by the energy-dependent broadening. 
By convolving the NBD and the probability for an inelastic collision over the impact parameter, the 
n-particle multiplicity distributions can be obtained. 

2. Energy Dependent Growth of Nucleon and Charged Hadron px Distributions 


In hadron-hadron collisions, the transverse momentum distributions for charged hadrons at lead¬ 
ing order can be expressed as [12] 
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where Nc = 3, xi ^2 = (PT/v^)®xp(±y), and Cp = (A^c “ l)/(2vr^A'c)- The running coupling con¬ 
stant 9in[Q2/A|^^] ; 0-5} with Aqcd = 0.2 GeV and = max{(kx + pr)^)/4, (kx - 

px)^)/4}. b and r^ are the impact factor and the transverse position of the gluon, respectively. 

In Eq.(l), (j) is the unintegrated gluon distribution (UGD). In the KLN model, cj) can be taken 
as [6] 
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where ^ is 


a normalization factor, and the saturation scale [4-5] 


( 3 ) 

where Qq = 2 GeV^, A = 0.288. In the KLR-AdS/GFT model, (/> can be obtained by a Fourier 
transform from the dipole-nucleus amplitude given in Ref. [8] 
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where the corresponding saturation scale 
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The parameter pm, Mq and Aq in Eq.(5), as given in Ref.[8], can be obtained by a fit to HERA data. 

In order to give an accurate theoretical analysis, the UGD should be considered as a function of 
the transverse position distribution of the nucleon. Thus, the gluon saturation momentum is always 
written as [4] 


<5s,p(A)(^>b) = Ql{x){ 
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where the nuclear thickness function 


^p(A)(b) = J ci2:pp(A)(b,, 


and Tp(A),o = 1-53 fm [4]. For the nuclear density distribution, we use the Woods-Saxon form [13] 


PAih,z) = 
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where r = Vb^ + and the measured values for Pb are po = 0.1612, ii = 6.62 fm, a = 0.545 fm. 
For the proton’s density distribution, the Gaussian form is used 
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Because of gluon saturation, the inelastic nucleon-nucleon cross section should grow as y/s in¬ 
creases. This will result in a broadening of the nucleon’s density distribution in position space. 
Therefore, the Gaussian width B should be written as a function of y/s [4] 
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with iTin(\/s) = 52, 60, 70.45, 72, 76.3 mb at y/s = 0.9, 2.36, 7, 8.8, 14 TeV, respectively [9,14-15]. 

In this paper, we also assume the gluon saturation scale should have a small dependence on y/s 
through the 3-dimensional rms radius of the proton [16] 


«b,A)(vs) = c5,p(A,(v;s)(^)'''t 


( 10 ) 


where S = 0.8 and the 3-dimensional rms radius rrms = /< r'^ > = /3B. 


3. Negative Binomial Distribution in the Glasma Flux Tube Picture 


Negative binomial multiplicity distribution is an interesting and important feature observed in 
multiplicity distribution of charged hadron production in high energy collisions. In the Glasma flux 
tube framework [11], the negative binomial distribution can be derived as 
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where the parameter k and mean multiplicity n are all considered as functions of the impact param¬ 
eter. The parameter k in the saturation approach is defined to be 


Hh) = C 


27r 


( 12 ) 


where C is a dimensionless parameter and S± is the overlap area of the two hadrons [11]. For a given 
impact parameter b, 

<3s,p'S'±(b) = J d2s_L(5^_p(s_L,b), (13) 

where (5s,p in the overlap area of the two hadrons is chosen to be (5s,p(s±, b) = min{(5s,p(s_L), (5s,p(s_L — 
b)}. The mean multiplicity h can be obtained by integrating Eq.(l) over 
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where the pre-factor Cm is proportional to cJm. 

In order to compute the probability distribution as a function of multiplicity, one should convolve 
the NBD at a given impact parameter (Eq.(ll)) with the probability for an inelastic collision (^^^^) 
at the same impact parameter. Thus, the probability distribution for producing n particle can be 
given by 

^ D . 

(15) 


P(n) = J d2b^pNBD(-(b),A;(b)). 


where the probability distribution can be given in impact parameter eikonal models [17] 

dPinel _ 1 6Xp( *^gg^pp) 
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with the energy dependent quantity Ugg = AttXB [4]. The overlap function for two protons at a given 
impact parameter can be expressed as 


Tp 
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4. Results and Discussion 

Charged hadron transverse momentum distributions in p-|-p collisions at different collision en¬ 
ergies are shown in Eig.l. Eig.l (a) and (b) are the results with the KLN and the KLR-AdS/CET 
model, respectively. The pT distributions are averaged over the ry range from -2.4 to 2.4. The rapidity 
(y) in Eq.(l) can be changed into the pseudo-rapidity (ry) using the transformation 

1 J cosh^ry -|- rriQ/p^ + sinhp 

y(ry) = -In , -, (18) 

J cosh^ry -|- — sinhp 

where mo(~ Aqcd) is the rest mass of particle. The experimental data come from CMS and ALICE 
[1-2]. The parameter ^(= 0.51) in Eq.(2) can be obtained by a analysis with the experimental 
data [18-19]. In Eig.l (a), it is shown that the theoretical results with the KLN model are in good 
agreement with the experimental data. For the KLR-AdS/CFT model, as shown in Fig.l (b), the 
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theoretical results fit to the experimental only at \/i = 7 TeV while pT is small. The reason is that 
the KLR-AdS/CFT model is an effective model at small Bjorken-rc {x < 10“^). At large px or small 
\/S) xi, 2 {= PtIV^ ■ exp(±y)) is out of the valid range of this model. In Fig.2, the predictive results 
with the KLN model for charged hadron px distributions in p+Pb (solid curves) and p+p (dashed 
curves) collisions are also given. Fig.2 (a) and (b) are the results at y/s =5.02 TeV and 8.8 TeV, 
respectively. It is shown that the results for p+Pb collisions are closer to the values of p+p collisions 
at small px than those at large px. 

Fig.3 (a) shows the the probability distribution for an inelastic collision at b. The curves are the 
results at y/s = 0.9 TeV (solid curve), 2.36 TeV (dashed curve), 7 TeV (dash-dotted curve), and 14 
TeV (dotted curve). The results show that the probability distributions 27r6dPinei/d^bx at certain 
y/s have a sharply peaked distribution and the impact parameter b of the peak grows with increasing 
collision energy y/s. Fig.3 (b) shows the width parameter k versus impact parameter b. The figure 
captions are the same as those in Fig.3 (a). It is shown that the width parameter k becomes larger as 
y/s increases. This is because the number of flux tubes Q'gSx in k depends on the energy-dependent 
broadening thickness function Tp. 

The results for charged hadron multiplicity distributions in the range of |r/| < 0.5 and |? 7 | < 1 
are shown in Fig.4. The figures are the results at y/s = 0.9 TeV (a), 2.36 TeV (b), 7 TeV (c), and 
14 TeV (d). The solid curves are the results with the Glasma flux tube approach, and the dashed 
curves are the results without considering the nucleon’s transverse position distribution. The data 
come from ALICE [20]. The parameter = 0.05 is extracted from a fit to data at y/s = 2.36 TeV 
and used for all other energies. It is shown that the theoretical results with the Glasma flux tube 
approach are in good agreement with the experimental data at all collision energies. For the method 
without considering the transverse distribution, the theoretical results fit well to the data at ^/s = 
2.36 TeV but show deviations for other collisions at highest multiplicities. 

In summary, charged hadron px distributions in p+p collisions at various collision energies are 
studied with the CGC approach. By considering an energy-dependent broadening of the nucleon’s 
density distribution in position space, the theoretical results of the KLN model are in good agreement 
with the experimental data from CMS and ALICE. Then, the probability distribution for producing 
n particle at different pseudo-rapidity ranges are also studied in the picture of Glasma flux tube, 
and the results fit well to the data. The predictive results for px and multiplicity distributions in 
p+p and p+Pb collisions will be examined by the forthcoming experiment at the LHC. 
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Figure Captions 


Figurel: Charged hadron px distributions for p+p collisions averaged over the range of \r]\ < 2.4 
with the KLN model (a) and the KLR-AdS/CFT model (b). The experimental data come from 
CMS and ALICE [1-2], 

Figure2: The predictive results for charged hadron px distributions at y/s =5.02 TeV (a) and 
8.8 TeV (b). The sold and dashed curves are the results for p+Pb and p+p collisions, respectively. 

FigureS: The probability distribution for inelastic collision (a) and the width parameter k (b) as 
a function of impact parameter b. 

Figured: Multiplicity distributions of charged hadrons in the range of |r/| < 0.5 and |r/| < 1. The 
figures are the results at y/s = 0.9 TeV (a) , 2.36 TeV (b), 7 TeV (c) and 14 TeV (d). The solid 
curves are the results with the Glasma flux tube approach, and the dashed curves are the results 
without considering nucleon’s position distribution. The data come from ALICE [20]. 
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